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Abstract. We overview applications exponential asymptotics and analyzable 
function theory to difference equations, in defining an analog of the Painleve 
property for them and we sketch the conclusions with respect to the solvability 
properties of first order autonomous ones. It turns out that if the Painleve 
property is present the equations are explicitly solvable and in the contrary 
case, under further assumptions, the integrals of motion develop singularity 
barriers. We apply the method to the logistic map a' n +i = ax n (l — x n ) where 
it turns out that the only cases with the Painleve property are a = —2, 0, 2 and 
4 for which explicit solutions indeed exist; in the opposite case an associated 
conjugation map develops singularity barriers. 



1. Introduction 

We present an outline of new methods for determining the position and the type of 
singularities of a certain kind of solutions of difference equations [17] and use this 
information to perform Painleve analysis on them. The approach relies on advances 
in exponential asymptotics and the theory of analyzable functions [10, 11, 23, 13, 
14, 15, 16, 18, 20, 21, 25, 26, 34, 45]. The main concepts are discussed first. 

Analyzable functions. Introduced by Jean Ecalle, these are mostly analytic func- 
tions which at singular points are completely described by transseries, much in the 
same way as analytic functions are represented at regular points by convergent se- 
ries. In contrast with analytic functions (which are not closed under division) and 
with meromorphic functions (which fail to be stable under integration and compo- 
sition) analyzable ones are conjectured to be closed under all operations whence 
the grand picture of this theory, that all functions of natural origin are analyzable. 
In particular, solutions of many classes of differential and difference equations have 
been shown to be analyzable. 

Transseries. Also introduced by J. Ecalle, transseries represent the "ultimate" 
generalization of Taylor series. Transseries are formal asymptotic combinations of 
power series, exponentials and logarithms and contain a wealth of information not 
only about local but also about global behavior of functions. One of the simplest 
nontrivial examples of a transseries is 



Key words and phrases. Borel summation; Exponential asymptotics; Singularity analysis; 
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x In the sense stemming from Stokes original papers and the one favored in exponential asymp- 
totics literature, Stokes lines are those where a small exponential is purely real; on an antistokes 
line the exponential becomes purely oscillatory . 
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What distinguishes the expression (1.1) from a usual asymptotic expansion is the 
presence beyond all orders of the principal series of an exponentially small term 
which cannot be captured by classical Poincare asymptoticity. More general trans- 
series arising in generic ordinary differential and difference equations are doubly 
infinite series whose terms are powers of x multiplied by exponentials (sec e.g. 
(2.3)). These transseries can be determined algorithmically and usually diverge 
factorially, but can be shown to be Borel summable in a suitable sense in some sec- 
tor. In this sector the function to which the transseries sums has good analyticity 
properties and on the edges of the sector typically singularities appear. 

The correspondence between the formal universe and actual functions. Envisaged 
by Ecalle to be a "total" isomorphism implemented by generalized Borel summa- 
tion, the correspondence has been established rigorously in a number of problems 
including ODEs, difference equations and PDEs. 

2. Difference equations: the isolated movable singularity property 



The problem of intcgrability has been and is a subject of substantial research. 
In the context of differential equations there exist numerous effective criteria of 
integrability, among which a crucial role is played by the Painleve test (for extensive 
references see e.g. [33]). 

An analog of the Painleve property for difference equations turns out to be more 
delicate to define. 

In the context of solvability, various methods have been proposed by Joshi [38], 
Ablowitz et al. [1], Ramani and Grammaticos (see references in [36]), Conte and 
Musette [12]. See also [1] for a comparative discussion of the various approaches 
in the literature. None of these is a proper extension of the Painleve test. One 
difficulty resides in continuing the solutions x(n) of a difference equation, which are 
defined on a discrete set, to the complex plane of the independent variable n in a 
natural and effective fashion. The embedding of x(n) must be done in such a way 
that properties are preserved. It is important for the effectiveness of the analysis 
that this embedding x(n)\Zx(z) is natural, constructive and unique under proper 
conditions. 

It is of course crucial that we are given infinitely many values of the function 
x(n); since the accumulation point of n is infinity, the behavior of x at oo is key, and 
then the question boils down to when infinitely many values determine function, 
and in which class. 

To illustrate a point we start with a rather trivial example. Assume that x(n) 
is expressed as a convergent power series in powers of 1/n 
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We would then naturally define x(n) □ x(z) by 
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x(z) = ^c k z k 

k=l 

for large enough z. Uniqueness is ensured by the analyticity at oo of x(z). Still by 
analyticity, C preserves all properties. 

We cannot rely merely on analytic continuation since there are very few classes 
of equations where solutions are given by convergent power series. However, as 
discovered by Ecalle and proved in detail very recently in a general setting by 
Braaksma [11], a wide class of arbitrary order difference equations admit formal 
solutions as Borel summable transseries. The class considered by Braaksma is of 
the form 

(2.1) x„ + i = G(x„, n) = + ^Aj x„ + g(n, x„) 

where G analytic at oo in n and at in z„, under genericity assumptions [11]. In 
particular a nonresonance condition is imposed 

(2.2) [i m = k • mod 2iri 

with k g N™ iff k = e m . Transseries solutions for these equations have many 
similarities to transseries solutions of differential equations. With some mi, < 

mi < n, 

(2.3) x(n) = J2 C k e- k -' XIl n k - a y k (n) 

In (2.3) /x = (fi\, ...,fi mi ) and a = (ai, ...,a mi ) depend only on the recurrence (they 
are the eigenvalues of A, A, respectively), C is a free parameter and y k are formal 
series in negative integer powers of n, independent on C. The number mi is chosen 
so that all the exponentials in (2.3) tend to zero in the chosen sector. 

Braaksma showed that yk(^) are Borel summable uniformly in k. Let Yk = 
'C _1 yk(«)- Then Y^(p) are analytic in a neighborhood of R + (in fact in a larger 
sector). Defining 

/>oo 

(2.4) y k = / e- n PY k (p)dp 

Jo 

we have uniform estimates |yk| < A k and thus the series 

(2.5) x(n)= C k e- k ^"n ka yk(n) 

kGN™i 

is classically convergent for large enough n. Braaksma showed that x(n) is an actual 
solution of (2.1). It is natural to replace n by z in (2.4) and define: 

(2.6) x(z)= C k e- k -^z k - a y k (z) 

kGN m i 

If z and all constants are real and /Uj < 0, the functions (2.6) are special cases 
of Ecalle's analyzable functions. As explained before we are allowing for z and 
constants to be complex, under restriction 3?(k • fiz) > 0. 
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It is crucial that the values of y(n) for all large enough n uniquely determine 
the expansion. In [17] it is shown that under suitable conditions, two distinct 
analyzable functions cannot agree on a set of points accumulating at infinity. Below 
is a simplified version of a theorem in [17]. 

Theorem 1 ([17]). Assume 

(Z • n = mod 2m with ZeZ p )«Z = 

and x(z) = X)keN"*i C k e _k ' tz 2; k ' a y k (z). If x(n) = for all large enough n e N, 
then x(z) is identically zero. 

Analyzable functions behave in most respects as analytic functions. Among the 
common properties, particularly important is the uniqueness of the extension from 
sets with accumulation points, implying the principle of permanence of relations. 
Under the assumptions [11] that ensure Borel summability, we have the following. 

Definition 2. // 

x(n) = J2 C k e- k ^"n ka y k (n) 

keN m i 

then we call 

x(z) = J2 C k e- k -^ k ' a y k (z) 

keN m i 

the analyzable embedding of x(n) from N to a sector in C. 

Having now a suitable procedure of analytic continuation, we can define the isolated 
movable singularity property, an extension of the Painleve property to difference 
equations, in a natural fashion. In analogy with the case of differential equations 
we require that all solutions are free from "bad" movable singularities. 

Definition 3. A difference equation has the IMSP if all movable singularities of 
all its solutions are isolated. 

Notes, (i) We use the common convention that isolated singularity exclude 
branch points, clusters of poles and barriers of singularities. 

(ii) To determine the singularities, transasymptotic matching methods intro- 
duced for differential equations in [18], can be extended with little changes to dif- 
ference equations. 

3. Classification of some difference equations with IMSP. 

Solvability. 

We look at autonomous difference equations of the form x n+ \ = G(x n ) where G 
is meromorphic and has attracting fixed points. A more general analysis is given 
in [17]. We then write 

(3.1) x n+ i = G(x n ) := ax n + F(x n ) 

and restrict for simplicity to the case F(0) = F'(0) = and < \a\ < 1. There is 
a one-parameter family of solutions presented as simple transseries convergent for 
large enough n, of the form 
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(3.2) x n =x n {C)=Y,z nkX ™C k D k 

fe=i 

for given values of Dk, independent of C. The analyzable embedding of x, cf. 
Definition 2, reads 

oc 

(3.3) x{z) = x(z;C) ^^e zk ^ a C k F) k 

fc=i 

which is analytic for large enough z. To look for the IMSP, we find the properties 
of x(z) beyond the domain of convergence of (3.3), and find the singular points of 
x{z). 

Note. Because equation (3.1) is nonlinear, although (3.3) has one continuous 
parameter, there may be more solutions. This issue is addressed in ([17]). 

3.1. Embedding versus properties of the conjugation map. By the Poincare 
conjugation theorem applied to x n +i = G(x n ) there exists a unique map </> with 
the properties 

(3.4) 0(0) = 0, <j>'(0) = 1 and analytic at 

and such that x n = 4>(X n ) implies X n+ \ = aX n . The map is a conjugation map 
of £ ra +i with its linearization X n+ \. We have X n — a n Xo — Ca n . 
We obtain an extension of x from N to C through 

(3.5) x n = 4>{Ca n ) C x(z) := (f>{Ca z ) 
Then the conjugation map satisfies the equation 

4>{az) = G(z) = a<j>(z) + F((f>(z)) 

As expected by the uniqueness of the embedding in Definition 2 we have the 
following. 

Lemma 4. (i) For equations of the type (3.1) under the given assumptions, the 
embeddings through analyzability and via the conjugation map, cf. (3.5) agree. 

(ii) The solutions x(z;C) have only isolated movable singularities iff (j) has only 
isolated singularities. 

Corollary A necessary condition for (3.1) to have the IMSP is that the conjugation 
map 4>, around every stable fixed point of G and of those of its iterates which are 
meromorphic ( of the form G^" 1 ' — G o • • • o G m times) extends analytically into the 
complex plane, except for isolated singularities. 

3.2. Autonomous equations with the IMSP. In ([17]) we classify the equations 
of the form (3.1) with respect to the IMSP. In a way analog to the case of ODEs 
of first order, only Riccati equations have this property. On the other hand, these 
equations are explicitly solvable, and thus autonomous equations of the first order 
which have the IMSP can be solved in closed form. 
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Figure 1. Julia sets for G = ax(l — x) for a = 0.5, a = 0.9 
and a = 1 respectively. The set K p is the interior of these curves 
(graphs generated using C and MapleV) 



Theorem 5. The equation (3.1) under the given assumptions fails to have the 
IMSP unless, for some c e C, 

(3.6) *) = T^I 

In case (3.6 we have: l/x n = a~ n (C — c/(a — 1)) + c/(a — 1) 

It is also very interesting to see that in the case of failure of the IMSP the analytic 
properties of the solutions preclude the existence of nice constants of motion. This 
is discussed in the next section. 

3.3. Case of failure of IMSP. We illustrate this situation when G polynomial. 
The surprising conclusion is that for polynomial G without the IMSP, constants of 
motion (defined as functions C(x,n) constant along trajectories) develop barriers 
of singularities along some fractal closed curves dJC pi see below. In a neighborhood 
of the origin, the function <j> is invertible, and thus for small x we derive from 
x n = cf)(a n C) that C = a~ n Q(x n ) where Q — . 

In Fig. 1. we depict the Julia set dJC p of a simple map. In that case, in the 
compact set bounded by JC P consists in the initial conditions for which the solution 
of the iteration converges to zero. The Julia set is a closed curve of nontrivial fractal 
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dimension. For a comprehensive discussion of Julia sets and iterations of rational 
maps see [5]. 

Theorem 6. Assume G is a polynomial map with an attracting fixed point at the 
origin. Denote by JC P ( cf. Fig. 1 ) the maximal connected component of the origin 
in the Fatou set of G. 

Then the domain of analyticity of Q is K p , and dJC p is a barrier of singularities 

ofQ- 

The proofs of these results can be found in [17]. The logistic map discussed in 
relative detail in the next section represents a very simple illustration of some of 
the relevant phenomena. 

4. Analysis of the logistic map at the superstable fixed point infinity 

We show that the equation x n+i = ax n (l — x n ) has the IMSP iff a € {—2, 0, 2, 4}. 
The case a = needs no analysis. Otherwise, taking y = l/x we get 

(4.1) Vn+l = 



a{y n - 1) 

For small y , the leading order form of equation (4.1) is y n +\ = — a ^ 1 Vn whose 
solution is — y\ a~ 2 _1 . It is then convenient to seek solutions of (4.1) in the 
form y n = F{y 2 a~ 2 ) whence the initial condition implies F(0) = 0, F'(0) = —a. 
Denoting y 2 a~ 2 = z, the functional relation satisfied by F is 

(4.2) F(z 2 ) = - Jgg i} (F(0) = 0, F'(0) = -a) 

Lemma 7. (i) There exists a unique analytic function F in the neighborhood of 
the origin satisfying (4-2) and such that F(0) — and F'(0) = —a. This F has 
only isolated singularities in C if and only if a G {—2,2,4}. In the latter case, the 
equations (4-2) and (4-1) can be solved explicitly (see (4-5), (4-8) and (4-9)). 

If a £ {—2,2,4} then the unit disk is a barrier of singularities of F. 

(ii) If a £ {—2,2,4} then dJC p is a barrier of singularities of the constant of 
motion Q ~ F^ 1 (cf. Theorem 6). 

4.1. Proof of Lemma 7. 

4.1.1. Analyticity at zero. We write 



(4.3) F{z) = ^F(z 2 ) - ^a 2 F(z 2 ) 2 - AaF{z 2 ) 

(with the choice of branch consistent with F(0) — 0, F'(0) — —a), and take F(z) — 
—az + h(z). This leads to the equation for h 

(4.4) h(z) = az+^(h(z 2 )-az 2 )-^Aa 2 z 2 - Aah(z 2 ) + a 2 (h(z 2 ) - az 2 ) 2 = N{h) 

It is straightforward to show that, for small e, Af is contractive in the ball of radius 
\a\ 2 in the space of functions of the form h{z) = z 2 u(z), where u is analytic in 
the disk D = {z : ||z|| < e} with the norm \\h\\ = sup D \u\. The corresponding 
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F is analytic for small z and is a conjugation map between (4.1) and its small-?/ 
approximation. 

4.1.2. Now the question is to determine the singularities of F in the complex 
plane. Inside the unit disk we can inductively continue F analytically through 
(4.3) as follows. If F is analytic in a disk of radius r < 1 then (4.3) provides the 
analytic continuation in the disk of radius y/r if F(z 2 ) is not zero or 4/a. Because F 
is assumed to be analytic in the disk of radius r, F(z 2 ) cannot vanish; otherwise by 
(4.2) it would vanish infinitely often in a neighborhood of zero, which is inconsistent 
with its analyticity at z — and with F'(0) = —a. But the equation F(z 2 ) = 4/a 
does have solutions in the unit disk if \a\ > 5. Indeed, with \a\ > 5 it is immediate 
to show F(z 2 ) — 4/a implies F(z 2n ) := e n — > as n — > oo. Since F is, by the 
condition F'(0) = —a, bijective near z = 0, the equation F(z) = e n has a (unique) 
solution, z = zq for sufficiently large n. Now, if z\ is such that z\ = zq we have 
F(z ) = 4/a. But then F has a square root branch point at z\ since F'{z) cannot 
vanish anywhere inside the unit disk, otherwise again by (4.2), F'{z) would vanish 
infinitely often near the origin. 

Thus for \a\ > 5 there is no F with only isolated singularities. In fact, as is 
seen in the last paragraph of the proof, in this case the unit disk is a barrier of 
singularities of F. 

We now consider the case where \a\ < 5. 

4.1.3. We claim that unless a = —2, the point z — 1 is a singular point of F. 
Indeed, a Taylor series expansion F = J2T=o c k( z —^) k gives c = or c = a/ (a— 1). 
It is straightforward to see that c = implies c k = for all k which contradicts 
F'(0) = -a. 

Therefore Co = a /(a — I) in which case direct calculation shows that, unless 
a=— 2(2 fe — 1) for some k € Z, all for k > 1 are zero which is not possible since 
F(0) = 0. If k > 1 then \a\ > 5. 

Therefore, k = 1 whence a = —2. With this value of a, (4.2) has the explicit 
solution 

(4.5) F(z) = - 2 - 

z 2 + z + 1 

It remains to look at the cases when z — 1 is a singular point of F. 

4.1.4. If z = 1 is a mcromorphic point of F, then we obtain from (4.2): 

(4.6) km ^ = 1 

v ' *-i F(z) 

If F(z) = Y^k=- P c k( z ~~ l) fc i s tne Laurent series of F at z = 1, then (4.6) implies 

(4.7) a = 2 p 

and, since \a\ < 5, we must have a € {2,4}. 
For a — 2 (4.2) has the explicit solution 



(4.8) 
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For a = 4 the solution of (4.2) is 



(4.9) 



F{z) = 




4.1.5. If F is not meromorphic at z — 1 then F cannot be continued analytically 
beyond the unit circle: 

Assume first that F is analytic in the open unit disk D\. By (4.2) we have 



where R n is a rational function. Thus if z 2 " — 1 and F is analytic at zo then it is 
meromorphic at z = 1, while the set of points such that Zq = 1 is dense on the 
unit circle. 

The last case to analyze is R < 1 , where R is the (maximal) radius of analyticity 
of F. Equation (4.3) shows that on the disk of radius R the singularities of F are 
square root branch points. On the other hand, if z e D\ is an algebraic branch 
point and z\ — zq then z\ is also an algebraic branch point as can be seen by 
continuing (4.2) around a small circle near the branch point z 2 = zq. It is then 
easy to see that if R < 1 the branch points accumulate densely towards the unit 
circle. □ 

(ii) The proof follows the lines of ?? so we only sketch the main steps. We note 
that by (i) F is invertible near the origin, so Q is analytic near the origin, where it 
satisfies the relation Q 2 (z) = Q(z 2 /(az — 1). 
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